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Abstract 

The undeformed limit of the dilute two-color braid-monoid algebra gives a natural basis for 
the description of spin-1/2 ladder models, and allows different Baxterization Ansatze. Based 
on this observation we find the complete classification of exactly soluble generalized isotropic 
spin-1/2 ladder models. 

1 Introduction 

Exactly soluble models afford considerable insight into the complicated behavior of frustrated 
quasi-one-dimensional spin systems. Here we investigate the most general spin ladder which 
includes isotropic Heisenberg leg, rung and diagonal interactions, as well as biquadratic leg- 
leg, diagonal-diagonal and rung-rung spin exchange terms. The ladder system represents 
the minimal model admiting both multiple-spin exchange interactions and non-trivial exact 
solutions. 

The spin models on a ladder are usually characterized in terms of their symmetry, i.e., 
the invariance of the Hamiltonian with respect to transformations in ordinary space and in 
spin space. In Ref. the sit (4) and su(3|l)- invariant integrable Hamiltonians have been 
constructed on the basis of the solution of the Yang-Baxter equation (i?-matrix), given in 
terms of the corresponding permutation operator. Multiparametric models related to these 
symmetries have been considered in |2] and [3] while an extension to n-leg ladder models 
associated with the su(2 n ), o(2 n ) and sp(2 n ) algebras has been obtained in jlj and 0. The 
i?-matrix of the generalized ladder model 4, Q satisfy the Yang-Baxter equation with a mul- 
tiplicative composition law. Physical properties of some exactly soluble ladder systems, such 
as magnetisation plateaus and thermodynamics have been presented in [Q, jHj, |10j . 

From the general point of view, many known examples of solutions of the Yang-Baxter 
equation are based on algebraic structures, such as braid algebra, Temperley-Lieb and Hecke 
algebras, Birman-Wenzl-Murakami algebra. Therefore it is natural to search for the algebraic 
construction suitable for describing SU (2)-invariant spin ladders. 

We study exactly soluble models within the space of coupling parameters. The general 
ladder Hamiltonian is based on a plaquette algebra, by which is meant an algebra of operators 
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defined on the four-spin plaquette units of the ladder. This algebra satisfies quasi-local com- 
mutation relations, meaning that operators defined on next-to-nearest (and further-) neighbor 
plaquettes commute. Using the commutation relations of this algebra we find exact expressions 
for different i?-matrices which give solutions to the Yang-Baxter equation. These solutions 
correspond in turn to soluble (and SU (2)-invariant) spin-ladder Hamiltonians, which may be 
either gapless or gapped. 

The plaquette algebra provides a representation of the undeformed limit of a dilute, two- 
color braid-monoid algebra. The dilute, two-color braid-monoid, or two-color Birman-Wen- 
zl-Murakami algebra was introduced and investigated in a series of papers by Grimm and 
coworkers jllj-|15j. In particular situations, when the generators of the undeformed dilute 
algebra are combined into full-braid and monoid generators, the algebra is equivalent to the 
Brauer |16j or braid- monoid algebra , which is the undeformed limit of the Birman-Wenzl- 
Murakami algebra. By taking all possible different combinations of generators of the algebra 
we construct the solutions of the Yang-Baxter equation, and thus the corresponding exactly 
soluble Hamiltonians. Because the algebra provides a complete basis for SU{2) invariant 
ladder we believe that the list of isotropic SC/(2)-invariant integrable ladder models (with 
short-range interactions) presented here is complete. 

In Section 2 we introduce the general spin- 1/2 ladder model and express all isotropic in- 
teractions in the plaquette-algebra basis. By seeking commutation relations of this algebra 
we observe that it provides a representation of the undeformed limit of a dilute two-color 
braid-monoid algebra. In Section 3 we use the commutation relations of the algebra to es- 
tablish a number of Baxterization Ansatze by considering different subsets of the generators 
of the algebra. We thus obtain all exactly soluble Hamiltonians in the parameter space of 
the generalized ladder Hamiltonian. In Section 4 we point out the relation between words, 
defined on the plaquette algebra and Matrix-Product ansatz states. In the concluding section 
we summarize the results and make contact with known integrable models. 

2 The ladder model and its algebraic structure 
2.1 Model and operator basis 

The Hamiltonian of the two-leg ladder under consideration is 

H = Jl(Sj • Sj + i + Tj • Tj+i) + JrSi • Tj + Jd(Si • T i+ i + T, • S i+ i) 

i 

+V LL (S l • S i+1 )(T* • Tt+i) + V RR (Si • T,)(S m • Tf+i) + V DD (S l • T m )(S m • T f ), (1) 

where we include leg, (Jl), rung, (Jr), and diagonal (Jd) Heisenberg interactions, as well as 
leg-leg (Vll), rung-rung, (Vrr), and diagonal-diagonal, (Vdd), four-spin interactions. Spin- 
1/2 operators on rung i are represented respectively by Sj and Tj, i = 1, . . . , L, and periodic 
boundary conditions are assumed, Sl + i = Si and Tl+i = Ti. Many solutions listed in the 
next Section are invariant with respect to the intertwinning transformation 

&2i <-> &2i T^2i <-> T?2i T2j+1 S2i+1, 

Jl <-> Jd Vll <-> V DD (2) 
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which interchanges S and T on every second rung and therefore we present one representative 
Hamiltonian for each pair related by this transformation. 
We introduce the orthonormal basis 

75" n> 

|2) = ^(IU) + UT» |3> = |U> (3) 

for each rung, and construct a set of corresponding projection operators which generate the 
SU(4) algebra (hereafter Greek indices denote 0, 1, 2, 3) 
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The spin operators Sj and Sj may be expressed in terms of these projection operators. The 
ladder Hamiltonian is thus equivalent to a four-state chain with only nearest-neighbor 
interactions. 

We introduce the basis of operators 
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where Latin indices a, b take the values 1, 2, 3 and we have used the plaquette notation 
Oi t i + i = Oi for all operators. 

This basis provides the full set of operators for isotropic interactions, i.e. any SU{2)- 
invariant term may be represented by linear combinations of the above generators. The 
multiplication table of these operators has the form 
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where in a product of two operators of the form o\ o\ , the operator o\ is taken from the 

(2) 

column on the left whereas the operator 0\ belongs to the top row. 

In the operator basis defined by Eq. Ip|). the Hamiltonian Q has the form 

H = J2 9l B i + 92di + 93(bi + b\) + g 4 (ai + o|) + g 5 p® + g 6 p 3 + c, (6) 



where 



Jl Jd . Vll , Vdd Jl Jd , Vll , Vdd 

Ql = Q2 = 

y 2 2 8 8 y 228 8 

— — _ in _i_ _ Yee — — — In — Ym _i_ ^ dd 

93 ~ ~2 2~ ~8 8" 94 ~ ~2 2 8~ ~8~ 

Vrr Jr Jr , Vll , Vdd , 3Vrr 

55 = ~ + T 9& = -^2 + — + — + — 

Jr Vll Vdd 3Vrr 



4 16 16 16 



2.2 Operator algebra 



(7) 



The operators defined in © generate a closed algebra whose commutation relations for oper- 
ators defined on the same plaquette are given in the multiplication table, while for operators 
on neighboring plaquettes one obtains 



e;ej+iej = e;p° +1 e i+ ieie i+ i = p°e i+ i 

BiB i+ iBi = B i+ iBiB i+ i eiB i+ \ei = eip® +l 

Biei + i = Bi + iei€i + i etB i+ i = eiei + \Bi 
Biei + \Bi = B i+ ieiB i+ i (8) 

for operators ej, called monoids and Bi, called braid operators, and 

b Wi+l = b i+l4 b l4+l a i = b i+l e i 

b \ b \ + l = CLi+l4 b \ b \+l a i = a i+l e i 

b\e i+ i = b i+l a\a i+ i b\e i+1 bi = b i+ ieib\ +l 

b\B i+1 bi = b i+1 Bibl +1 (9) 

bib i+l a\ = 4+iPi hcii+iei = b\ +l ai 

ha i+ ia\ = p}b\ +l bia i+ iBi = b\ +1 aiB i+ i 

bib i+1 Bi = B i+1 bib i+1 (10) 

4h+\bi = (k4 + i 4 a i+i b i = e i b l+i 

4h+i = ei4+ibi 4 a i+i = e i b \+i b \ ( n ) 
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OioJ+i = h+ibi 
a-A+i = a i+ ibi 
a-i&i+i = ciiBi + iBi = 6j + i6jOi+i 



a A+i b \ 
a M+i b i 



bi+w] 



aiB i+ iei = aip9 +l 
o,i&i+\Bi = a,iBi + i 

aie i+1 a\ = aiB i+ ia\ = pfp] +l 



(12) 



ejaj +1 = a\b i+1 bi = B i+1 B t a\ +1 e^e^al = a\p\ +1 



eib, 



(13) 



-BjaT +1 = Bj + ieioJ +1 Bia\ +1 b\ = B i+1 a[bi + x 
Bie i+ ia\ = B i+ xa\ BiB i+1 a\ = e i+1 a} 
BibUA = bt +1 b\B i+1 . (14) 

between operators of different kind. For \i — j\ > 1 all operators commute. 

This algebra provides a representation of the dilute two-color braid-monoid algebra intro- 
duced and investigated in ^J-JHl) ° r ) more explicitly, to its undeformed limit. It emerges 
in the investigation of dilute A-D-E models, and may be considered as a special case of the 
two-color generalization of the braid-monoid algebra ^7]) where the representation of one of 
the colors is trivial. This algebra possesses a graphical interpretation in terms of strands 
or strings, such that the relations of the algebra are equivalent to the deformation of the 
diagrams. 

The operators e« satisfy the relations of the Temperley-Lieb algebra ^B] (see the first line 
in Eq. ©)• The braid generators Bi, together with e«, satisfy the braid-monoid (or Brauer 
jlHj ) algebra (jSJ , a subalgebra of the dilute two-color braid- monoid algebra. Later we will also 
introduce full braid and full monoid generators which satisfy the 0(4)-related Brauer algebra 
relations. All of these algebraic relations permit the construction of different Baxterization 
patterns. 



2.3 Baxterization 

The Yang-Baxter equation plays a central role in the study of integrable models. By starting 
from the braid-group representation one may construct its solutions, which depend on the 
spectral parameter, and demand that in the limit, as the spectral parameter is taken to be 
zero, the given braid-group representation is recovered. The process of seeking this solution is 
called Baxterization. It was introduced by Jones in Ref. |19j . where it was found that the Bir- 
man-Wenzl-Murakami algebra provides the solution of Yang-Baxter equation. In the present 
considerations, the algebraic relations ([9H14|) permit the construction of different Baxterization 
patterns. 
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The solution of the Yang-Baxter equation 



Ri(u)R i+ i(u + v)Ri(y) = R i+ i(v)Ri(u + v)R i+ i(u), (15) 

where u is real, is sought in the form 

Ri(u) = a(u)(ai + oj) + 6(«)(6j + b\) + c(u)Bi + d(u)ei 

+e(u)p° + /(u)^ 1 +p?) + g(u)pl (16) 

Following the standard analysis we require the initial condition, Ri(0) = Ii, and the unitarity 
condition, Ri(u)Ri(—u) = k(u)k(—u)Ii, for some function k(u) which we demand to be either 
rational, trigonometric or hyperbolic (in this work). Following the algebraic Bethe ansatz 
technique (see e.g. Ref. [20 for a review), the corresponding soluble Hamiltonian is constructed 
as 

^ = E^u>. (it) 

i 

This equation amounts to a strong condition on the general Hamitonian © and thus restricts 
the values of coupling parameters, for which the model is integrable, to a few points and lines 
in parameter space. 



3 Exact solutions 

We study different particular cases by taking certain operator combinations. The set of "basic" 
generators {B{, ej, Gtj, a\, 6j, &J} form a basis for all interaction terms in Q except for rung and 
rung-rung interactions. The hermiticity of the Hamiltonian requires that operators (oj, a\) 
and (bi, b\) enter only as the sums dj + a\ and bi + b\. Simple combinatorial arguments give 
the total possible number of exactly soluble models. Four combinations may be constructed 
from single generators (i.e. Bi, ej, a% + a] and b{ + b\), there are six pair combinations, four 
triple combinations and one combination constructed from all four generators. All of these 
cases must be considered together with the projectors pf. The combination a% + a\ enters 
only together with the monoid e^, as otherwise the unitarity condition for the R matrix is 
not satisfied. This excludes one single, two pairs and one triple combination. Thus there are 
in total 11 possible combinations which may potentially yield exactly soluble models, but we 
will show that not all of these are relevant to the spin-1/2 ladder or possess a Baxterization 
Ansatz. For example, two triple combinations, {Bi,bi + b\,ei} and {Bi,cii + oj,ej}, appear 
not to be integrable. In these cases the integrability test for the local Hamiltonian ha+i, 

+ [hi,i+i, hi+ij+z]] = 0, (18) 

i 

is not satisfied, whereas for all the integrable cases investigated in the present work it is 
satisfied. 

For certain cases we will present particular discrete transformations which increases the 
number of soluble Hamitonians. 
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3.1 The combination pf} 

This combination gives the solution to the Yang-Baxter equation for which the corresponding 
i?-matrix has the form 

Ri{u) = Ii + uik + bl+pl+pf). (19) 

This solution has exactly the form of the spin-1/2 Heisenberg model. Indeed, the composite 
operator IT = h + b\ + p ® + pf satisfies the braid algebra relation Iljllj+illj = ni+illjllj+i 
and its square gives the identity operator. The Hamiltonian which follows from this i?-matrix 
has exactly the same eigenspectrum as the spin-1/2 Heisenberg model, but the degeneracy 
of this spectrum is different 22 . The corresponding Hamiltonian -f/i/2 = J2i^i commutes 
with the generators J2i-^i b f° r an y CL,b = 1,2,3. These generators form an su(3) subalgebra 
within the su(4) algebra generated by the X-operators. i?i/2 commutes with the operator 
Y^i Xf°, which represents the total number of rung singlets, and any multiple of this term may 
be added to the Hamiltonian without affecting the integrability. It corresponds to a coupling 
of the effective spin-1/2 Heisenberg model to an effective "magnetic field" h. It is well known 
that the spin-1/2 chain in a magnetic field develops an incommensurate critical phase for 
\h\ < 2 and has a massive phase for \h\ > 2. Returning to the original spin variables, we have 

H1/2 = ^ Sj • Sj.fi + Tj • Tj + i — Sj • Tj + i — Tj • Sj+i 

i 

+4[(S< • S <+ i)(Ti • T i+ i) + (Si • Ti)(S i+ i • T i+1 ) 

-(Si ■ T i+1 )(S m • T«)] + (Jr + 2)(S< • Ti), (20) 

where J/j = correspond to the Hamiltonian derived from the ^-matrix (|19j) . The Hamil- 
tonian commutes with ^^Xf , and remains therefore integrable for arbitrary J/j. For —4 < 
Jr < 4 the model is in the critical, incommensurate phase, while for Jr > 4 it has a gapped 
rung-singlet phase, and in the region Jr < —4 it has a rung-triplet ground state. It coincides 
with the model studied in Ref. [22] and more recently in Ref. 0. Although the i?-matrix 
of Ref. 6 is not the same as in (|19|). we believe that they are related by an appropriate 
transformation. 

In the critical region the model is described by a conformal field theory with central charge 
c = 1. However, because of the degeneracy this conformal field theory possesses additional zero 
modes, and therefore the conformal dimensions appearing in this model have the Coulomb-gas 
dimensions of the c = 1 theory, but with different degeneracies. 

3.2 The combination {bi, b\, B i: pf } 

This combination leads to the ^-matrix 

Ri{u) = h + u(±(6j + &{) + Bi+ pf ), (21) 

a solution which corresponds simply to local Hamiltonian operators proportional to the permu- 
tation operators (see Sec. 3.5 below) pf (for positive sign) and P* (for negative sign). These 
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Hamiltonians are invariant with respect to the full 5^7(4) group, and describe the soluble 
"spin-orbital" model [231 

-^Pll = ' + ^ ' r ^ i+1 + ^(^* ' Si+l)(T» • Tj + i). (22) 

j 

The model corresponding to the Pj X can be obtained from the (|22|) by applying the transfor- 
mation ((2j) . Both Hamiltonians commute with J2i Xf° , and therefore the term Jr J2i Sj • Tj 
may be added to these Hamiltonians without affecting the property of integrability. The 
model (1221) was studied in Ref. jTf for certain values of the rung-rung interaction parameter 
Vrr added to the Hamiltonian. For some regions in the parameter space of Jr and Vrr, 
these models demonstrate a variety of gapless behavior. In particular, the model (|22[) in the 
continuum limit corresponds to the SU(4) Wess-Zumino-Witten (WZW) model at level 1 of 
the Kac- Moody algebra [21], [IS]- The central charge is equal to 3. 

3.3 The combination b\, e^, p®} 

This yields the hyperbolic solution 

~ , . sinh(A — u) , n q, , i 9s sinh(u) . „ . , 

R i( u ) = sin V h(A) ' (Pi + Pi) + (Pi + Pi) + (ei + <K&i + &!)), (23) 

with cosh(A) = 3/2 and <r = ±1. 

The corresponding Hamiltonian is 

H = — Sj • Tj + i — Tj • Si+i + 4[(Sj • Sj + i)(Tj • Tj+i) 

i 

-(Si • T;)(S i+1 • T i+1 )} + J R (S l • Ti) (24) 

for <t = +1 and = —1, and the case of a = — 1 is obtained by the transformation @. We 
note that these Hamiltonians commute with the rung singlet number operator J2iXf°, and 
that the model is therefore integrable for any Jr. 

3.4 The combination {bi, ej, a\,pf} 

This combination also permits Baxterization. The solution is given by the Izergin-Korepin 
model, which is related to the dilute A 2 model |26j.|27j. 
6 _ sin(2A - u) sin(3A - u) sin(3A - u) x 2 sin(u) sin(3A - u) 3 

n i u — • fn\\ ■ in\\ Pi ' ■ /o\\ \Pi 'Pi) ' v 1 ' ■ /o\\ • fo\\ 'Pi 

sm(2A) sm(3A) sm(3A) sm(2A) sm(3A) 

sin(u) sin(3A — u) „ + , . sin(u) sin(A — u) sin(w) , + . , . 

sm(2A) sm(3A) sm(2A) sm(3A) sm(3A) 

However, the value of the Temperley-Lieb factor (or fugacity in corresponding statistical me- 
chanics models) n in the relation ef = nei is restricted here by the condition n = — 2cos(4A). 
In our case n = 3 (see multiplication table) and therefore A is necessary complex. The analytic 
continuation of the weights into the domain \n\ > 2 by the change A — > 7r/4 + iX and u — > iu 
makes these complex, and the corresponding Hamiltonian then has complex coefficients which 
make this operator combination not relevant for the spin-1/2 ladder model. 
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3.5 The full set of operators 



The full set of operators provides solutions related to the Brauer algebra l$\ , or undeformed 
braid-monoid algebra |17j . In order to present this solution we introduce full braid and monoid 
operators, whose relations with the generators © of the dilute algebra are given by 

E\ ~ ] = ei + ai + o} + Pi E\ +) =ei- ai - a\ + p? 
pf =B i + b i + bj+pf P* =B l -b l - b\ + pf (26) 

where PjL_i and P* i+ i are respectively the permutation operators corresponding to leg- leg 
and diagonal-diagonal bonds. 

The operators pj^li, p|i~+i, P\ i+1 and P t x i+ i are generators of the algebra 



Ef = ^Ei EiEi + iEi = Ei Ei + \EiEi + i = Pj +1 

Pf = 1 PiPi+lPi = Pi+lPiPi+l P%Ei = EiPi = Ei 

PiEi+iEi = Pi+iEi Ei + \EiPi + \ = Ei + \Pi PiEi + \Pi = Pj+iPjPj+i 

and 

P%Pj = PjPi EiEj = EjEi PiEj = EjPi 

for |i — j| > 1. (27) 

We have again used the plaquette notation P^i+i = P% and P^j+i = Ei, where Pi denotes 
either pj' or P* and Ei denotes either e\ or e\ + \ Note that there exist a number of 
relations between the operators which follow from the above. The operators Pj satisfy the 
braid algebra relations, while the operators Ei form a Temperley-Lieb algebra [T^]- Together 
these operators form a Brauer JH] (or braid monoid J2j) algebra. The operators E± are 
projectors on the plaquette singlet states (below). 

This algebra allows Baxterization J7],[2S], the solution for the case under consideration 
being given by the rational R matrix 

u 

Ri{u) = h + uPi - ^jEi, (28) 

which satisfies the initial and unitarity conditions with k(u) = 1 + u. The spin Hamiltonian 
which corresponds to this R matrix is 

H = J2(Pi ~ Ei). (29) 

i 

We thus obtain four different soluble Hamiltonians which correspond to the four combinations 
{PJ',PH}, {pJ',p| +) }, {P*,PH}, and {P- x ,pf +) }. Two of these are trivially soluble, 

because X)i(pJ' ~~ E^) is the Hamiltonian of two decoupled chains and J2i{Pt X ~ e\ ^) is the 
same pair of chains intertwined by the transformation l[2*]). Of the remaining two, one is given 
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by the combination {Pj ' , e\ '}, which in terms of spin operators is 

Hb = Sj ■ Sj + i + Tj • Tj+i + Sj • Tj + i + Tj • Sj + i 

i 

+4[(S 4 • S m )(T* • T i+ i) - (S, t • Ti+i)(Si+i • Ti)] (30) 

The last, {Pj X , E^" 1 ^}, is obtained from the previous Hamiltonian by the same intertwining 
transformation (J2J). 

The Brauer algebra (|2*T|) is related to the representations of the group D2 = 0(4), and 
therefore the corresponding Bethe Ansatz is also related to this algebra. In fact the corre- 
sponding operator C is the product of two ^-operators for 6- vertex models. The eigenvalues of 
the algebraic Bethe Ansatz are given by the product of two sets of eigenvalues corresponding 
to 6- vertex models and finally one obtains two decoupled Bethe Ansatz equations which 
correspond to the D2 = SU(2) x SU(2) algebra. The corresponding energy spectrum is the 
sum of the eigenvalues for each SU (2) component. Thus the Hamiltonian (|3U|) is equivalent to 
the Hamiltonian of two decoupled chains. This model is critical (no spin gap) with a central 
charge c = 2. 

The two projectors e\ ' and E^ are related by the Aj-operator transformation 

A° a -» -iX° a Xf -> iXf (31) 
for all sites i = 1, ..,N. This is a unitary transformation generated by the operator 

N 

U(ir/2)=eM-4j2(X?°)}- (32) 
z i=l 

It leaves P, invariant and transforms the Hamiltonian of two decoupled chains into Hb ■ More- 
over, because of the relation [C/(7r/2)] 4 = 1, this transformation is one of the generators of the 
Z4 group. 

The projectors may be decomposed into the product of two operators, namely 

4~ } = A\A t E\ +) = B\Bi, (33) 

where 

A i - A j A i+1+ A i A i+1 ~ A « A i+l _A i A i+1 

a i — A i A i+l _A i A i+1+ A i A i+1+ A i A i+1- K 6 ^) 

These operators obey the relations 

Aj = At Bf = Bi 

M = 4X°°X°° 1 B t B\ = 4Af A™ . (35) 



The transformation l)32|) yields A{ <-> B{ and e\ ■* <-> P 4 - + ^ . The unnormalized projectors may 
be expressed as 

Et } = iv^xv^i 4 +) = \4 +) )(4 +) i (36) 
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where 

|^H)=4|0) |^ + )>=SJ|0>, (37) 

and we have defined the vacuum state as a direct product of singlet states for all rungs. The 
states \ip\ ^) and |V>j- ) may be expressed in terms of local singlet states = 1, 2, 3, 4) 

on the bonds of each plaquette. Using the X-operator notation and denoting the plaquette 
vertices in clockwise order, one may write 

\4^) = -2(\s l2 )\s u )-\s 23 )\ SA1 )) \4^) = -2\s 23 )\ SA1 ). (38) 

These functions yield a plaquette-singlet state, i.e. the state with the total spin per plaquette 
(Si + S2 + S3 + S4) 2 equal to zero. There are exactly two possible plaquette-singlet states 
which correspond respectively to e\ ' and e\ . The Z4 transformation (|H2|) maps one state 
into the other. 

In terms of the original spin variables this transformation is given by [HO] 

S i = ^(S i + Ti)-S < xT i T^KSi + T^ + SiXT; (39) 

The transformation is canonical for the spin-1/2 ® spin-1/2 representation of SU(2) x SU{2) 
group and the values of the Casimir operators are the same, S| = T? = 3/4. In terms of the 
variables Sj,Tj the Hamiltonian Hq (|3U() indeed has the form of two decoupled spin chains. 

3.6 The full monoid generators e\ ' and e\ + ^ 

These combinations allow Baxterization according to the Ansatz for the Temperley-Lieb al- 
gebra, 

~ . . sinh(A — u) r sinh(u) _ . , . 

where cosh(A) = 2. Here E\ denotes either e\ ^ or e\ + ^ defined in Eq. (|2(ij) . The corresponding 
Hamiltonians are 

H = - Y,{Si • T m + T, • S i+ i - 4(S, • T m )(S m • T<)} (41) 

i 

for the full Temperley-Lieb generator e\ , and 

H = - ]T{S, ; • S i+ i + Ti • T l+1 - 4(S, • S i+1 )(T i • T m )} (42) 

i 

for generator e\ . We note that these Hamiltonian operators are related by a transfor- 
mation. In this particular case two Hamiltonians (|41[) and (|42j) may also be related by the 
transformation ([%]). The spectrum of both Temperley-Lieb soluble models is gapped and their 
ground-state energy and gap may in principle be calculated by a mapping onto the Potts 
model |3*T] . 
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3.7 Models related to the bilinear-biquadratic spin-1 chain 

In this subsection we project the space of states onto a three-state subspace by applying the 
projection operator p® on every site. This gives three soluble models which are related to 
the generalized spin-1 chain. The role of unit operator Ii is played by the operator p®. The 
realization of these models in the spin-1/2 ladder is discussed below. 

3.7.1 The combination {B^e^p®} 

This combination of the braids and monoids of the dilute algebra ® form the Brauer algebra 
1(2*71) in a restricted 3-state subspace with Temperley-Lieb (fugacity) factor n = 3 (in the 
relation ef = net), and leads to the solution 

ik{u) = if + uB t - ^y$H- (43) 

This .R-matrix solution corresponds to the Takhtajan-Babudjian spin-1 integrable model [32| . 
The corresponding Hamiltonian is 

3 

i 

-2[(S< • Si+iXT; • T i+1 ) + (Si • T <+ i)(S i+ i • TO]. 

(44) 

This model is described by a c = 3/2, SU(2)f : =2 WZW model in the continuum limit. 

3.7.2 The dilute braid generator Bi 

This generator is a permutation operator in the 3-state subspace and leads to the SU{3)- 
invariant Uimin-Lai- Sutherland [33] model, 

Ri(u) = 4 3) + uBi. (45) 

The spin Hamiltonian is 

i 

+2[(Sj • S i+ i)(Ti • T i+ i) + (Si • T i+1 )(S i+1 • Ti)]. (46) 
In the continuum limit it is described by the SU(3)k=i WZW model |3i] . 

3.7.3 The dilute monoid operator 

This operator generates the Temperley-Lieb algebra in the 3-state restricted subspace (see 
relations ©, (fl3|0 . and leads to the solution given by the usual Baxterization Ansatz for the 
Temperley-Lieb algebra, 

Ri(u) = ^Z^/f) + !^ e? , (47) 
smh A smh A 
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where cosh(A) = 3/2. The corresponding Hamiltonian 

H = - £ Si ■ S t+1 - T, • T m - S, • T m - T, • S m 
+4[(S< • S m )(Ti • Tt+i) + (Si • T i+1 )(S w • TO] 

(48) 

is related to the spin-1 chain model studied in Ref. |31| . 
3.7.4 Generalization 

The three solutions above correspond to soluble, spin-1, bilinear-biquadratic chains. Their 
relation to the ladder (with four state per rung) is the following: on rung i the effective spin- 
1 variable is formed by the generators X? , (a = 1,2,3). If there is a finite region on the 
ladder formed by rung triplets bounded by rung singlets, X^^Xf X l ^_ l ....X'jp +i X^ Sri+1 , one 
may consider this region as an effective spin-1 chain with open boundary conditions. The three 
ladder Hamiltonians of this subsection thus yield exactly soluble models for such configuration. 

An explicit relationship between the original spin-1/2 variables and effective spin-1 oper- 
ators can also be given by the composite-spin representation [3"Kj . This approach was used in 
Ref. [HI. 

Considering the ladder with alternating domains of rung triplets and rung singlets, one 
notices that the couplings Vrr and Jr are arbitrary, because the Hamiltonian commutes not 
only with J2i Xf° but also with Xf° for any i, and therefore the eigenstates of the correspond- 
ing ladder Hamiltonian in ^-dimensional state space are characterized by the value of the 
total spin at each rung of the ladder. One may study the dependence of the ground state 
energy as a function of the Jr, Vrr and the total number of rung-singlet bonds. As shown in 
Ref. [SHI depending on these parameters there exist several gapped phases in such a ladder. 

4 Ground state and relation with Matrix-Product 
Ansatz 

We mention further possible implementations of the dilute algebra used in this analysis. The 
Matrix-Product Ansatz was proven in Ref. jSZ] to be a valuable technique for the description 
of the gapped states of spin systems. It is based on a decomposition of the spin Hamiltonian 
into the sum of projection operators on plaquette states with fixed angular momentum. There 
are two plaquette states with total momentum equal to (j = 0), three triplet states (j = 1), 
and one quintuplet (j = 2). For particular combinations of coupling constants it is possible 
to find exact ground states. In the current formalism, states with j = are created by linear 
combinations of the generators pf and a\ acting on the product of rung-singlet states, states 
with j = 1 are created by the operators bi ± b\, and by the combinations B{ — and p 1 ± p 2 , 
and the state with j = 2 is given by 6(Bi + p\ — 4ej). These states have the form of the 
Matrix-Product Ansatz |57| . 
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A natural interpretation also exists for these states as "words" , the letters of which are 
elements of the dilute algebra defined by the operators in Eq. (pjh The relation between 
exact ground states and words of a (Temperley-Lieb) algebra was noticed first in Ref. |38| . It 
provides a direct algorithm to search for possible exact ground states: one first classifies all 
possible simple words of the type i"Q_iLjMj + i, where in this analysis the operators K,L,M 
may be any of ai, aj, 6j, &j, -Bj, and pf, and then constructs a product of these elementary 
words. If the Hamiltonian leaves this product invariant we have obtained an exact eigenstate. 
In general this procedure can be used for constructing variational wave functions. Examples 
of such states are given by the following words 

\1>rs) PipI-Pn-xPn and \^ FM ) — ► p\pI-Pn-iPn (49) 
which correspond to the rung-singlet and ferromagnetic ground states respectively, and 

N 

|^=o) = II lVQ=o>i, (50) 
1=1 

where the local state 

ICi=o>* = (1 + 3^1/2 (Pf~ x 4M (51) 

produces the variational wave function which interpolates between the rung-singlet state (for 
x = 0) and one of the plaquette-singlet states (for x = 1). 

5 Conclusion 

We have shown that the undeformed dilute two-color braid monoid algebra provides a natural 
tool for the description of integrable spin- 1/2 isotropic ladder models with isotropic nearest 
and next-nearest-neighbor interactions. The Baxterization of different subsets of generators 
of the algebra yields corresponding spin Hamiltonians which are non-trivial. Their general 
feature is the presence of multiple-spin exchange, which may in some cases lead to critical 
behavior. The "words" constructed as a sequence of the elements of the dilute algebra yield 
variational wave functions of the ground state. 

At this point it is worthwhile to compare the present work to related studies of integrable 
ladder models. In fact, some of the models obtained above are not new. For example the 
spin-orbital type model of Section 3.2 is known from the work on soluble SC/(4)-invariant 
spin ladders. The class of models related to the integrable bilinear-biquadratic spin-1 chain 
(Section 3.7) has also been constructed using the composite-spin representation [33], jSHJ. A 
spin Hamiltonian, similar to Eq. lj20j) has already appeared in Ref. jJJ, but the corresponding 
i?-matrix is different. We believe that the connection between these two different /^-matrices 
can be understood on the basis of Hecke-soluble models of Ref. Additional support of 

this argument is given by Ref. |5|. 

In contrast to the models of Sections 3.1, 3.2 and 3.7, which have appeared before in 
the literature, those of Sections 3.3, 3.5, and 3.6 appear to be new. Recently, Batchelor and 
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coworkers 5j have used operator algebras for constructing soluble n-leg spin ladders in terms of 
o(2 n ) and sp{2 n ) -invariant it!-matrices based on Temperley-Lieb and Birman-Wenzl-Murakami 
algebras. But their construction of algebra operators is different from ours. For the two-leg 
ladder their Temperley-Lieb generator Q has a fugacity factor —4 (Q 2 = — 4Q), while our 
generator E has a fugacity factor equal to +4 (see Sections 3.5 and 3.6 ). Correspondingly, the 
resulting spin Hamiltonians are different as well. We also notice that although the Hamiltonian 
(|3()|) may be obtained by a ^-transformation from the Hamiltonian of two decoupled spin- 1/2 
chains, the algebraic structure underlying its integrability was not known before. 

Applications of these exactly soluble models to real spin ladders may be made by the use 
of field-theoretical methods. The exact solutions which describe gapless points in parameter 
space correspond to certain conformal field theories, WZW models, in the continuum limit. 
The perturbation of these theories by relevant and marginal operators which correspond to 
different terms in a real spin Hamiltonian produce in general a renormalization-group flow away 
from the gapless behavior thus allowing the determination of the properties of the relevant 
phases in the vicinity of the points of second-order phase transitions. This program is realized 
in the context of real ladder systems with cyclic four spin-exchange interactions in Ref . jHS] • 
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